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NATIOKAL ADVISORY COMMITIEE FOR AEROKAUTICS 


RESEARCH MEMORAUDUM 


THEORETICAL INVESTIGATION OF LAMINAR HEAT TRANSFER ON 
YAWED INFINITE CYLINDERS IN SUPERSONIC FLOW AND 
A COMPARISON WITH EXPERIMENTAL DATA 
By Ivan E. Beckwith 

SUMMARY 


A theoretical method is presented for the calculation of heat trans- 
fer in the laminar boundary layer on yawed infinite cylinders in conipres- 
sible flow with a Prandtl number near unity. The method may be applied 
to a cylinder of arbitrary cross section and arbitrary chordwise wall- 
temperature distribution. General expressions are derived for the heat 
transfer and recovery factor in the vicinity of the stagnation line . The 
expression for the recovery factor and the chordwise velocity profile at 
the stagnation line for zero heat transfer are in good agreement with the 
results of other theoretical analyses . 

.i 

The method is appl^^pWifeeit^he calculation of local heat-transfer rates 
on a yawed circular cylinder at a stream Mach number of 6 . 9 . Comparison 
of the results of this eiRcfSlai^c^wit!^ some experimenteil data indicates 
good agreement for small yaw ah0.es but at the larger yaw angles of 60 ° 
or more the theoreticaU-jjjpredicted heat transfer is about 40 percent less 
than the measiired values. 

INTRODUCTION 


As flight speeds are increased, the problem of aerodynamic heating 
becomes more acute and requires that every possible means be investigated 
for alleviating the heat transfer at critical locations on missile or 
aircraft configurations. One such criticeil location is a sharp leadi n g 
edge or nose where the local heat transfer becomes very large because of 
the thin boundary layer which has just begun to form. Appreciable reduc- 
tions in heat transfer to the forward regions of wings and bodies may 
therefore be obtained by rounding or blunting the leading edge or nose. 
(Another advantage of a blunt entry section over a sharp section is that 
heat may be removed more effectively from the blunt section by heat 
conduction or internal cooling) . Recent theoretical and experimental 
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investigations (refs. 1 to 4) as well as early experimental resiolts 
(ref. 5) liave indicated that fiirther reductions in heat transfer to hlunt 
leading edges can he obtained hy yawing or sweeping the wing. 

Inspection of the 1 nar houndary-layer equations for a yawed 
infinite cylinder in an inconpressihle flow with a Prandtl nxanber of 
unity and negligible viscous heating shows thht both the chordwise flow 
and the temperattire distribution are not affected by yawing the cylinder 
(ref. 1), that is, the tenperature and chordwise velocity distributions 
are independent of the spanwise flow component. The solution for this 
case (ref. 2) then shows that the heat transfer decreases as the sq.uare 
root of the cosine of the yaw angle. This decrease in heat transfer is 
due to the increase in boundary-layer thickness which occurs at a given 
point on the cylinder as the yaw angle is increased. 

In a compressible flow the independence principle for the yawed, 
infinite cylinder no longer applies since shear stresses associated with 
the spanwise flow will affect the temperature distribution directly, as 
well as affecting the chordwise velocity profiles throvigh the resiilting 
density changes, (refs. 6 to 8). An analytic solution of the general 
problem of heat transfer to yawed infinite cylinders in compressible flow 
is at present vii*tually impossible so that recourse must be had to an 
approximate solution such as can be obtained from a suitable extension 
of the integral methods (refs. 9 to 12 ). The range of validity and 
accuracy of a solution of this type cannot be determined without com- 
parison with mathematically exact solutions or experimental data. As far 
as is known, the only exact solution available for the compressible flow 
problem is Crabtree's (ref. 8) in which th&isitoCferential equations for 
Prandtl nmber unity and zero heat transfer are solved in a linearized 
form. , • 

The purpose of the present paper is to present a theoretical inves- 
tigation of the compressible laminar boundary'' layer on yawed infinite 
cylinders and to compare the results with some e^erlmental data showing 
the effect of yaw on the heat transfer to a ■simulated leading edge 
(ref. 5)- The theoretical analysis is based on an integral method wherein 
a linear viscosity temperature relation is used and the boundary layer is 
assumed to have one thickness. The latter assxmption simplifies the cal- 
cvilations considerably since the simultaneous solution of only two differ- 
ential eqiiatlons is required to obtain the local heat-transfer rates. A 
more general method would Involve the simultaneous solution of three 
differential equations; one for the thermal boundary-layer thickness, and 
two more for the boundary-layer .thicknesses corresponding to the chord- 
wise and spanwise components of the velocity bovtndary layer. 
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SYMBOLS 

"bovindary-layer coordinate parallel to surface, normal to 
, cylinder axis 

Boundary-layer coordinate parallel to surface and cylinder axis 

boundary-layer coordinate normal to cylinder surface 

time 

mass density 
velocity in x-direction 
velocity in y-direction 
velocity in z -direction 
pressure 

dynamic viscosity 
kinematic viscosity 
static enthalpy, CpT 

thennal conductivity 

static temperature , 

stagnation temperature, °R 

recovery temperature for local q:^ = 0 

stagnation enthalpy (eq. 6) 

Prandtl number , 

specific heat at constant pressure 

constant in equation of state for gases 
Dorodnitzyn variab^ ■*'10 ) ^ 
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P*- 

A 


’>*-z 


local stagnation density 

boundEiry-layer tMckness in t] coordinate system 


t^ = H - = H - CpOV 


X t*’ ^ ~ 

t = ~y— = £ — j thermal -profile parameter 


V 


Hi - 




S Sutherland constant 

A chordwise velocity-profile parameter (eq^. (I5)) 

Bi, B2j - coefficients in thermal profile (eq. (18)) 

0 momentum thickness in tj coordinate system (eq. (22)) 

6 * displacement thickness in t] coordinate system (eq. (25)) 

A* thermal -profile parameter (eq. (26)) 

9 energy- loss-profile parameter (eq. (27)) 

Uq, o^, Og • . . constants in expression for 9 (eq. (28)) 

Mq_ local chordwise Mach number external to boundary layer 




local heat-transfer rate at wall (negative for heat 
transfer to wall) 


L 


characteristic dimension of body 


velocity component normal to cylinder axis on stagnation 
streamline ahead of cylinder 

resultant velocity ahead of bow shock (flight velocity 
of cylinder) 


K> 

X = x/l 


resultant Mach number eihead of bow shock 
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h. 



yaw angle of cylinder (measured from a line normal to 
resultant stream direction) 


function of A and. 



(eq. (45)) 


chordwlse con^jonent of Mach number on stagnation streamline 
ahead of cy lin der 


pressiare coefficient 


diameter of circulax cylinder 

ratio of specific heats, 1.400 used throughout 

local heat -transfer coefficient 

average heat-transfer coefficient 



reference Reynolds number, 


Ps 


Z, 0 computing parameters (eqs. (A5) and (a4)) 

Subscripts : 

1 local conditions external to boundary layer unless 

otherwise noted 

w conditions at wall 

00 conditions ahead of bow shock 

B conditions just behind normal bow shock 

1q conditions external to boundary layer at stagnation line 


• ■■ 
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THEORY 

Derivation of the Integral Eg.uations 


The houndary-layer equations for a three-dimensional flow axe 
obtained by application of the conventional boundary-layer assumptions 
to the general equations governing the motion of a conipressible, vis- 
cous, heat-conducting gas (refs- 7 15) • Two of these assumptions 

which may be emphasized for the present application are (l) the boundary- 
]^yer thickness is small compared with some characteristic dimension of 
the body, and (2) the pressiore gradient normal to the surface may be 
neglected provided that the boundary-layer thickness is small compared 
with the local radius of curvature of the sijrface. Thus, for example, 
the boundary- layer equations are expected to be valid everywhere on a 
circular cylinder except in a separated region provided that the stream 
Reynolds number based on diameter (or radius of, curvature of the 
leading edge) is large. 

The boundary-layer equations for steady flow on an infinite yawed 

cylinder are obtained from the three-dimensional equations with -^ = 0 

^t 

and -^ = 0. The boundary-layer coordinate system used is defined by the 
following sketch; 



The eq-uation of continuity is then 
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and tbe eq.nations of motion in the x-, y-, and z-directions are 

written eis 


pu^ + = -k + 

Sx ^z ^ Sz V 5z/ 


( 2 ) 


pu^ + 
Sx 


^Z Sz V bz' 


( 5 ) 


0 = 


^z 


w 


where z is normal to the surface. The energy equation is written as 


Sh Sh 

pu — + pw — 

ax az 


. + I.(,M + ^ 

ax hz\hzj 


^au 

^az 


l\2 + 

J \^zj 


( 5 ) 


Then, for convenience, a stagnation enthalpy is defined as 


H = h + A 

2 2 


( 6 ) 


whereupon equations (2), (5), and ( 5 ) may he combined to yield 


pu^ + p„M = A/pM'\ + (1 _ i\ Afp^'\ 

ax az az \ azy / az \ az/ ■ 


( 7 ) 


for constant Cp and a. The chordwlse pressure gradient impressed on 

the botindary layer is obtained from eqmtion (2) evaluated outside the 
boundary layer, that is, outside the region where viscous effects are 
important, as 
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Pl^- 


dU]^ 

dx 


dx 


(8) 


which Is • Bemoul 1 1. * s eqixation. The eq.uation of state for a perfect gas 


p = pET 


( 9 ) 


is \ised throTighout. ^ 

Equations ( 2 ) and (7) are integrated across the houndaory layer (as 
in ref. 11) aind hy using equations (l), ( 4 ), (8), and (9), together 
with the Dorodnitzyn VEirlahle 


Tl = 




( 10 ) 


the results may he ejtpressed as 


_d^ 

dx 


A 




+ 




dT|* + 


A 




^1 dx 


^yPw / 8 _u\ 


(11) 


and 


ni 

— A -^{1 - t)dii* 
dx Jo ui 



\^1 dx 


t^* dx J 



^(1 - t)dTI^ = 


^H^Pw /8t \ 

(tpi^iA\8W^ 


(12) 
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for the houDf^-axy conditions at t]* = 0 of u = 0^ v = 0, w = 0, and 

¥ = 0 and at tj* = 1, u = ui, v = vi, t = 1, and — ^ = 0 where A 

Sti* 

is a universal houndaiy- layer thickness. Note that these eq^tlons Eire 
of exactly the same form as the corresponding equations of reference 11 
for two-dimensional flow. The spanwise component of the flow, however, 
is involved directly in the density-ratio integral of equation (ll) sind 
in the stagnation enthalpy Integrals of equation (12) as will he shown 
later . 


Velocity and Enthalpy Profiles 

An approximate solution to equations (ll) and (12) becomes possible 
after specifying the velocity and enthalpy profiles. The chordwlse 
velocity profile will be described by the conventional fourth-degree 

polynomial in q* = ^. A fourth-degree polynomial is then used for the 

A 

spanwise velocity profile in order to satisfy the requirement that both 
profiles must be of the san^ shape in the limiting case of a yawed flat 
plate. The stagnation enthalpy profile will be expressed as a fifth- 
degree polynomial which has resulted in acceptable engineering accuracy 
for the heat transfer to an unyawed cylinder when T^ < Tg (ref. ll) . 

The coefficients of a l l profiles are determined from specified conditions 
at the edge of the boundary layer (tj = A) and at the wa ll (q = O) so that 
the original equations of motion (eqs. (l) to (T))^ after being transformed 
to the new variable t], are satisfied when evaluated at the wall . 

The profile coefficients can be simplified if a viscosity-temperature 
relation of the form 



(15) 


is assumed where the quantity — is a function only of x with 

^w 

determined from the given values of T^j- by any acceptable method such as 
Sutherland's viscosity-temperature equation. In the following derivation, 
equation (15) is actually used in the vicinity of the wall only,, and, as 
will be seen, the final results in the form of a heat-transfer coefficient, 
for example, are automatically expressed in terms of and -where 

may also be evaluated from T^^ according to the best available infor- 
mation. In the practical application of these results or for the purpose 

‘'(7. ^ y'T I' •, 
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of coraparison with experimental data, however, it is usually desirable 
to use a reference viscosity and conductivity such as and k^,. 

Consequently, the problem arises as to the proper evalmtion of the 
ratios According to the above considerations, and 

since aind . never enter the basic derivation, the viscosity 

ratio should be obtained from the Sutherland equation which would be 
i^ritten as- 


VV S 


'Dir conductivity ratio should also be determined from equally acctirate 
infomnation. 

The chordwise velocity profile may now be written as (ref. 11) 


JL = (14) 

5 2 -'6 6 

which satisfies the conditions at Tj* = 1 that 



The relation for A is obta±ned from equation (2) (transformed to i] 
and evaluated at t] = 0) as 
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A 




du- 


dx 


1 a2 


(15) 


after also iising equations ( 4 ), (8), (9), and (15) • 

Application of the same procedure to equation (5) shows that^ at 

Ti = 0, 


82 

8it*2 



= 0 


(16) 


Hence , the required spanwise profile can he written as 


i = 2(,n*) - + (n*)** 


which at tj* = 1 gives 



(17) 


Coraparison of equation (17) with equation (l 4 ) shows that^ for a flat 
plate (?v = 0), the two profiles become identical as would be expected. 
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The stagnation enthalpy profile is written as 


= t = (i8) 


where the boundary conditions at n* = 1 are taken as 


t = 1 


^ = 0 


with the result that 


B 4 = 5 - tel - 5B2 - 2 B 5 


( 19 ) 


B^ = —4 + 5®! + 2Bg + Bj 


The coefficient B 2 is derived, from eq.uation (j) which is transformed to 

the X, T) coordinates and evaluated at the wall (u = 0, w = 0). The 
result can be expressed as 


B. 




( 20 ) 


after uRi ng eq^uations (4), (9)^ ( 15 ) definitions of H and 

t. The first derivatives of the velocity profiles are obtained from 
equations (l4) and ( 17 ). 
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The third thermal profile coefficient is derived from the 

first derivative with respect to of equation ( 7 ) (expressed in the 
X, Ti coordinate system) evaluated at the wall. In order to obtain 

/ 

the final expression for Bz, as given below, note that — ) =0 

. \^n/q=o 

and (^1 ’ = 0 from equations (l) and (lO), respectively. Also 
equations (I 5 ) and (20) have been used to supply expressions for 

. The result may be written as 

T]=0 


and the quantity 




B = ^ i 12A + 

ti* 36 

\ ' / T1=0 


where derivatives of the velocity profiles from eqioations (l4) and ( 17 ) 
have been lased where necessary. The coefficient B]_ in equation ( 18 ) 

will be determined from the solution of the integral equations. 


dx 


- 3(1 - CT)U, 


( 21 ) 


Expression of Profile Integrals in Terms of X and Bj_ 


The various profile Integrals in equations (U) and (12) can now be 
evaluated in terms of X and Bj_ which are the two basic parameters to 

be determined from the simiiltaneous solution of these two equations. 


Thus, 
in the q 


from equation (lA) and the definition of a momentxmi thickness 
coordinate system, there is obtained 


9 

A 





A _ 

5 ihhj 


(22) 


Similarly, a displacenent thickness in the x, q plane (analogous to 
the displacement thickness in incompressible flow) is expressed as 
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M. = 

A 



^6 - A 
120 


(25) 


Then from the definition of H and. 
density integral may he expressed, as 


and the relation 


P Ti" 


the 



Tl \Ti* 


iW + e ^ 

J A 2CpTi A 



2CpTijA 



2CpTi 




(24) 


hy means of eqviationh (22) and (25). The last Integral on the right in 
eqiiation (24) becomes hy virtue of equation (17) 



(25) 


The integral expression ^ 
tions (18) and (19) as 


in equation 


(24) is evaluated from equa- 


^ = 
A 


(1 - t)dT]^^ " I ■ ^ 




(26) 


where equations (20) and .(21) are used, to express B2 and. as 

functions of A. The integral term in eqimtion (12) is evalmted in 
the same way and yields 
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£ = 
A 




The values of the a's depend only on the type of thermal profile used 
and are listed helow for the fifth-degree profile used here. 


Oq = 0.58175 

= 0.04484 

ot^ = -0.14605 

= -0.01250 

Og = -0.05517 

og = -0.00557 

= -0.01589 

= -0.00079 


(28) 


Eq,uatlons (ll) and (l2) may now be written in a more convenient form 
by introducing equations (22), (25), (24), (25), (26), and (27); thus, 


M + J- 

dx U]_ 


duj^ 

dx 

- 1 


2 + 7 “ ^ ^ 

“ 

M + 265 A 
0 650 0 2CpTi 


2W ^ 

2 / 0 

^^^2w 12 + ?v 

UjIpi^^ 6 


(29) 


and 

^ i _ _2 p_ 

dx, \^uj^ dx t^* dx^ o(pifuiA- 


(50) 
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where 


M 


L 


]ItR\ 


The external spanwlse velocity may be expressed as a function of 

temperature and the local chordwise Mach number M]_ as 


2CpTl 




(51) 


In principle^ equations (29) an.<i (50) can now be solved for ~K and 
by a suitable numerical procedTxre for any given distribution of 

and Tj^. Before the confutation can be started, however, the initial 
values and first derivatives at x = 0 of A and Bq_ must be known. 
(See appendix. ) 

General Expression for the Local Heat Transfer 
The local convective heat-transfer rate is defined as 


<3w = 



Introducing the transfonnatlon to q* and using equations ( 15 ) and (l8) 
and the definition of t* resiilt in 





V n/A 


(32) 


A heat-transfer parameter which depends only on the yaw angle for given 
stream conditions and wall temperature is then obtained from equation ( 32 ) 
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where for a symmetrical cylinder with detached how wave 


Ps %3 
cos 3 


P u 

*^00 CO 


This latter relation would be only approximate for an asymmetrical 
cylinder . - 


Application to Stagnation Type of Flow 

Initial values of A and In the vicinity of the stagnation 

lihe which occurs on a blunt leading edge for both the yawed and vmyawed 
cylinder^ u^^ approaches 0 as x approaches 0 . Then, in order to avoid 

infinite valiies of — and ^ at x = 0, the tvro following conditions, 

dx dx 

as obtained from equations (29) and (50), must be satisfied at x = 0: 


2i + 8^!: + ]_W + 265/V 

A A Ti^\ti* ;a 630\Tiq 



12 + A = 0 

Tio 6A 


( 5 ^) 


and 


£ 1 ^ 5 l = 

A CT Tiq ?v 


( 55 ) 


where equation (15) bas been used to supply a relation between A^, 

dx 

and Tv. The expressions for — , and ^ are Introduced ijato • 

AAA A 

equations (^k) and (35) from equations (22), (23), (26), and (27) with 
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B 2 and (from eq.mtions (20) and (21)) evaluated at x = 0. The 

parameter B^^ may then he eliminated between equations (5^) and (55) 

and the following quadratic for results; 




+ <crA 


20 + ^ 


A - 4 


2A 


Ti* 


(»i + + 5(00 + 


A A 


v% 


- 1 


2 ^ + JL(i . a) 
650 15 V / 


“■w 


‘io 




fTi* 


- 1 


<y) 

630 15 ' 


’I 


^A(l - a)^a2 + - 1^ = 0 


( 56 ) 


where for U]_ = 0 


Tn* 


T 


1 - 

‘lO 


= 1 + 


"1 _ 


1 + 


7 - 






1 + 


7 - 1 


cOS^ 3 


(57) 


since 


Vi^ = = 7RTJ1„2 gjji2 p 
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and 






1 + 


Tben^ for given values of M and p, equation (56) is solved for 


i 

Tin 


as a fimction of A. Typical results for tMs calculation are shown in 
figure 1 for various valijes of corresponding values of B]_ 

may he found from equation (55) s-s 



where equations (20), (21), and (27) evaluated at x = 0 have been used 

to obtain 2. 

A 

Recovery temperat\ire and heat-transfer coefficient at x = 0 .- For 
the condition of zero heat transfer, 0^ = 0 and the wall temperature 

eqiials the recovery temperature Tg. Combining equations ( 58 ) and (52) 
and using these conditions give an esgoresslon for Tg which is 


Tg = Tj^* + l4-(l - a) 


02 + ^ , 

°0 + 


T 


To, 


(59) 


Then defining a recovery factor r in terms of the local teiig)erature 
outside the boundary layer results in 
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r = 


Te -Tio 


Ti* - 




= X + 4(1 - (i) 




“s 

oq + |a4 


m 


■which, gives for -val-ues of a near unity r =» /o~ since equation (4o) is 
aljnost independent of (Actually, for zero heat transfer -fche -value 

of Tv must satisfy eqmtions (39) anii (3^) 'with = Tg.) This -value 

of the recovery factor is in good agreement -with the res-ults of an exact 
analysis hy Schuh. (ref. l4). 

A heat-transfer parameter using the conventional definition for the 
heat-transfer coefficient 


l^c = 


Sr 




m 


may now he formed hy substituting equations ( 38 ) and ( 39 ) Into equa- 
tion ( 33 ) so that 


hgL 


^ , 


PcoUcoL 




Tiq uTv 



(^2) 


where all quantities are e-valua-bed at x = 0. It ’may he no-fced that this 
relation is- expected -bo aijply in -the -vicinity of the stagnation line for 

an isothermal cylinder since ( — 1 = 0 and 

<^Jx=0 


in the appendix and lim 

X — >0 


G 

^ 

^ WL = o ' 


&Bi\ 

=0 as shown 

V'^/x=0 


The function 


F Tv, 


Oq + 




^In/ \ 




( 43 ) 
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appeaxing in eq.uation (te) has 1)6611 coniput6d with a = O .7 for a larga 
rang6 of th6 rasults ara plottad agalost in fig- 

11 T 6 2. This figure shows that ^ regardless of the values of and p, 

the function F depends almost entirely on Tjj^^T]_^ where 


n- 

^1q ^1* 1 + 2LzJm^2 cos2 P 


m 


from eq.uation (37). Thus, a numerical expression for F as a function 

T„ 

of — and p could he used to give a universal correlation for 

Tj_* 

the heat transfer along the stagnation line. Such an expression for 
a = 0.70 is 


FA, 




= 0.8k ^ 



- 0.5 


0.4 


(^5) 


with a maximum error pf about 4 percent for 2.5 > > 0.10. (See 

fig. 2.) The resulting correlation parameter from equations (42), (44) 
and ( 45 ) is then 


h(»L 


^ I 


P UcoL 


/cos 




0.84 

1ti* 


1 + 




- 0.5 


1 + ^ ~ COS^ py 


- 0.4 


(^6) 


which is a function of the yaw angle p, the stream Mach number 
the wall-to-sta®iation- temperature ratio T^^T]^*. The quantity 

depends on the chordwise conponent of the Mach number just behind the 
shock according to the isentropic relation 
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It shoTold be noted at this point that by analogy wltb the results of 
reference 11 equation (^6) is expected to underestimate the heat transfer 
coefficient by as much as 10 to 30 percent when > Tg (see fig. 5, 

ref. 11). 

Application to a 'circular cylinder .- The velocity-gradient parameter 

in equation (46) is evaluated from a modified Newtonian type pressure 
coefficient (ref. 15 ) which is written for a circular cylinder as 

= g-i u- = oos^i^ 

1 2 1 2 V T) 

— pu^ — pu'=^ 

2 00 00 2 OO 

where D is the cylinder diameter ♦ This eq^uation modifies the usual 
expression for the Newtonian pressure coefficient (ref. l6) so that the 
correct pressure coefficient is obtained at the stagnation line. An 
analogous equation may be easily derived for any other cross-sectional 
shape. Differentiat in g the equation for a circular cylinder twice with 
respect to x and evaluating at x = 0 give 



d%^ 

dx2 




m 


Bemoxilli's equation (8) requires that 


lim 

X— ^0 



d2 


■Pi 


<3x^ 


1/2 


Jx=0 




hence, from equation (48), 







NACA EM L55F09 


23 


or Introducing the chordwlse velocity behind the shock gives 




T- 


lo 


-|l/2 


Y 


(50) 


This relation has been plotted against the crossflow Mach num- 
ber Mw = Moo cos p in figure 3 where it is compared with some experi- 
^*00 

mental results on cylinders . The experimental values shown are the 
results of some previously unpublished data obtained in the Gas Dynamic 
Branch of the Langley Laboratory. The length of the vertical line 
through the data points corresponds to the maximum experimental error. 


Chordwlse veldcity profile at the stagnation line for zero heat 
transfer . - As an indication of the accuracy of the present method for the 
condition of zero heat transfer, the chordwlse velocity profile at x = 0 
is compared with the resiolts of Crabtree's linearized solution of the 
same problem (ref. 8). The chordwlse profile given by eq^uation (l4) is 
transformed back to the z -plane by means of the relation 


I 1 ^^1 % ^ 

Aviq ds Tp* 



Pn* 

-i- dTyif- 
P 


(51) 


which is the Inverse of equation (10 ) combined with the expression for A 

from equation (15) evaluated for T„. = Tn* and — ^ = —1— corresponding 

% % 


to the conditions of zero heat transfer, a = 1, and the viscosity- 
temperature relation of reference 8. The density ratio in the above 
equation, for these conditions, is written as 
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from the definitions of H and t and equations ( 4 ), ( 9 ), ( 1 ?)^ and 
(3T)- Inserting this density ratio in the first equation and integrating 
gives the final required relation between a and q*-. The value of A 
required in equations (l4) and ( 5 I) was obtained from equation ( 36 ) for 
these same conditions and is 7.249- 


The results of this calculation for a spanwise Mach number of imity 
are shown in figure 4 where a/Uq is plotted against the distance param- 


eter 




The reasonable agreement between the present method 


and Crabtree's more exact analysis indicates that the compressibility 
effects, in this particular case, eire predicted satisfactorily by the 
approximate integral solution. 


COMPARISON OF THEORY WITH EXPERIMENTAL DATA 


The experimental, data of reference 3 sxe representative of an average 
heat-transfer rate per \mit area or, alternately, the total heat-transfer 
rate for the forward half of a circular cylinder. A dimensionless param- 

SpD 

— was found to correlate the data for eill yaw 


eter of the form 


angles as shown in reference 3 * Ihe theoretical analysis applied to the 
local heat-transfer coefficient at the stagnation line yields for this 
same parameter the expression 






/ 


P 


1 

cos P — I 
Hg dXl 


'Ps 


0.84 


1 + 




.-0.5 


1 + 


7 - 


^2 


- 0.4 


COS^ 


(52) . 
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[X 

from equation ( 46 ) witli (Cp and a being assumed constant) 

and X = — . This has been plotted in the lower paxt of figure 5 for 
D 

T 

= 6.9, Tt* = 1200° R and —2. = 0.2, 0-5, 0.8, corresponding approxi- 

T’ * 

■^1 

mately to the experimental test conditions ^except for = 0 . 2 ^. Com- 
parison of the/ theory with the data indicates that the meas\ired effect of 
yaw on the average or total heat transfer on a half -cylinder is consider- 
ably less than the theoretically predicted effect on the local heat trans- 
fer in the vicinity of the stagnation line. Thus, in order to obtain a 
more coniplete comparison of the theory with experiment, it is necessary 
to calcTilate from the theoiy the local heat transfer over the entire half- 
cylinder. This calcxilation was carried out (according to a general proce- 
dure outlined in the appendix) for two angles of yaw, 3=0° and 3 = 60°, 

with M„ = 6.9, dT^/dx = 0 , i = 0.5 and 0 . 2 , and Tp* = 1200 ° R. 


dM*i 

The values of M-i and — ^ required in the calculation were obtained 

dx 

from the faired curves shown in figure 6. These curves are based on the 


/< 3 M. 




experimental data Included in this figure as well as values of 1 _ 1 

V^/x=0 

computed from the data shown in figure 5 * The data points shown in fig- 
ure 6 were computed directly from measured pressure data according to the 
Isentropic relation 


Ml = 


7-1 



- 1 


which, from inviscid theory for yawed infinite cylinders, would depend 
only on the val\ie of Mw . Consequently, since some of the data from 

reference I5 show considerable scatter, some additional data (previously 
unpublished) at a comparable have been Included. Then largely on 

the basis of these data for Mm = 4.08 and the results of figure (3), the 

*^00 

erratic data for O.5 > x > 0 could be overlooked and the ciirves shown 
result . 

j •' ' 
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The resiilts of the calciilation are shown in figure 7 where the ratio 
I ^* 3 ^ j ^ plotted against x for the two yaw angles and temperature 


/x =0 

ratios considered. Note that the variation of 


‘^/(‘^)x=0 


with yaw angle 


and temperature ratio is not very large and, for practical purposes, coifLd 
be considered a function of x only. Also Included in figure 7 "tbe 

QirD 

initial values of the basic theoretical parameter ^ 




from eq_uation (33)* The present method fails for A > 12 (see appendix) 
which occurred in all cases before x = ^ as indicated in figure 7. 


However, sufficient information was obtained to permit a reasonable 
extrapolation to ^ ^ which was used as the basis for the following 

calculations . 


The average heat-transfer coefficient h^, was obtained by integra- 


tion over X of the local cop5)uted values of and then division by 


(’V - tJ 


where T^ was taken as the experimental value according to 


the upper part of figure 5 
then used to conrpute the valiies of the parameter 


These average heat-transfer coefficients were 

h^ 

kco 




00 

,UcJ> 


which are 


shown in figure 8 . (Accurate values of m^/Mw ^/Ko were used in 

computing this parameter) . The experimental data in the form of this 
same parameter (ref. 3) sjre also plotted in figure 8 . Comparison of the 
theoretical values with the experimental data reveals that the predicted 
heat transfer is about percent less than the experimental results for 
a yaw angle of 60 ° whereas at a yaw angle of OP the agreement is good. 

The reasons for the discrepancy at large yaw angles are not yet clearj 
however, at least part of the difference may be due to the configuration 
of the test model used in reference 3 - The finite length and end shape 
would both tend to make the test cylinder approximate something between 
an infinite yawed cylinder emd a pointed body of revolution at angle of 
attack so that the measured heat transfer could possibly depend on the 
distance from the end as well as the yaw angle. Actually, in view of the 
approximating ass\mptions Involved in the theory (mainly, the universal 
boundary-layer thickness and the linear viscosity relation) and the possi- 
bility of certain indeterminate factors in the experiment (such as end 
conditions and spanwlse heat conduction), the agreement between theory 
and experiment is probably as good as can be expected. 


On the basis of the results of the present method at the stagnation 
line (eq. 40 ) and Schuh's theoretical resxalts (ref. l 4 ), the value of a 
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I 



"local" recovery factor would be 


r 


Te - ^1 

Ti^ - Ti 




(55) 


This assumption cannot be compared directly with the experimental valiies 
of reference 3 which represent an average recovery factor or tg[uilibrium 
/temperature for the condition of zero total heat input to the entire half- 
cylinder. However ^ an average equilibrium temperature may be defined and 
evaluated from equation (53) with the result that 



= /ff + (1 - /a)- 


1 + 




cos 






rl 


Jo ^lo 


die 


( 5 ^) 


where eqimtion (37) IIS'S been vised and depends only on M]_ from 

figure 6. This average which is a fvinction -of the yaw angle only 

for a given M<jo, is compared with the experimental values in the vpper 
part of figure 5 . Although exact agreement between theory and experiment 
is not expected since the local heat transfer is not everywhere zero in 
the experiment, the use of equation (55) is confirmed in this case. 


CONCLUSICNS 


An integral method has been developed for the calculation of the 
compressible l ami na.r boundary layer on yawed infinite cylinders of arbi- 
trary shape and with arbitrary wall temperature distribution. In general, 
the method requires the numerical solution of two simultaneous ordinary 
first-order differential equations . At the stagnation line of a blunt 
cylinder the solution reduces to a set of algebraic equations which result 
in a general correlation equation for the heat transfer in the vicinity of 
the stagnation line. An expression is also derived for the recovery fac- 
tor which is in good eigreement with an exact solution and some experimen- 
tal data. The chordwise velocity profile at the stagnation line as com- 
puted from the present method for zero heat transfer and a Prandtl number 
of unity is also in good agreement with a more exact solution. 
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The method Is applied to the calculation of local heat-transfer 
rates on the front part of a circiilar cylinder at a stream Mach numher 
of 6.9 and angles of yaw of 0° £ind 6o°. The results of this calculation 
indicate that the variation in local heat transfer around the surface 
of the cylinder has .very little functional dependence on yaw angle or 
wall temperature i Comparison of the calculations with some experimental 
data, which shows the effect of yaw on the average heat transfer to a 
simulated leading^ edge , gives good agreement at small yaw angles hut at 
a yaw angle of 6o° the theoretical heat transfer is about 40 percent less 
than the measured value. This discrepancy may he caused partly hy inde- 
terminate factors in the data as well as the simplifying assumptions 
involved in the theory. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 7 , 1955- 
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APPENDIX 

CALCULATION PROCEDURE 
General Method 


Equations (29) and (30) are written in the form 


^ + 2Z 
dx 


2 + 


r - 


iM 2 A 2U:^ 2W aifi . 
1 V 2 1 / 0 Tq^ 




26^ 

63 o| 


•Tl V 2 1 y 


► — ^ = 2-^^^ 12_^ 0 (AIJ 

' dx p^* u^p^*L “L 6 


and 


M + 2(il^ ^ S ^w 

dx ^Iut dx ti* dx / Cf rv* U,D,** 




m 


(A2) 


where the parameters Z and 0 


have heen introduced so as to avoid the 


iise of the second velocity derivative 


d^ui 

dx^ 


which is difficult to obtain 


accurately .from experimental data. The qijantlty Rg^^ is the Reynolds 

number at some reference point, for example. Just behind the bow shock on 
a blunt cylinder. The computing parameters are defined as 
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and 


JjL 

dui 
dx 

where equation ( 15 ) has been used to provide the relation between A and 
A. For convenience, all quantities needed for the solution of these 
eqioations are listed below: 





0 = 

A 315 




(22) 


^ = 56 - A 
A 120 


( 25 ) 


^ It'D X -p X - p 

T ■ 3 ■ 5^ ■ ^ 5 


(26) 


£ 

A 


- °0 + 6 °^ ■+ (°1 + ^)®l + ^“2 + 6°^)^ (°5 6^)®5 


Bg = (1 - 0 ) 


Ti^ - OV 




- ^ 


ig-pr -M 


(A3) 


®5 


= 2L^ 2_a_ 122 l^ L - 1 + ^ 

o X m Jf. m A 3 ,, 




Ti* - 6 


3ui ^ 


dX; 


(a6) 
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■ 


^ Ti Tiq 



1 + 

7 - 1 p p 

-_-i^ cos2 ^ 


(AT) 


where eqiiation (5?) has heen iised in the last expression. Equations (A5) 
and (a 6) are obtained from equations (20) and (21) with equation ( 5 I) used 
where necessary. The simultaneous solution of equations (Al), (A2), (A5), 
and (Ah) for Z, "K, 0, and can now be carried out by any suitable 

graphical or numerical procedure. The values of u-j^ ^or Tj^, 

Moo, and p must be given. A typical numerical procedure is applied to 
the same type of problem in reference 12. 


Initial Values . 


Blunt leat^ing edge .- In order to start the integration the values at 
X = 0 of Z, Tv, 0 , Bq^, and their first derivatives must be known. Eor 
a blunt cylinder the initial values of Tv and B]_ are obtained from 

equation ( 36 ) (or fig. 1 ) and equation ( 58 )* corresponding values of 

Z and 0 may then be obtained from equations (A3) and (a 4) . For the 

condition that ~ and ^ are finite at x = 0 on a blmt cylinder, 
dx dx 

it can be shown from eq^uations ( 29 ) and ( 30 ) that (^1 = 0 and 


dBq\ 

<^x/x=0 

eqiiatlons 


= 0 provided that 
(A3) and (a 4), it 


V dx^ 


= 0 and 


'x=0 

follows that ( — 
dx 



= 0 


x=0 


x=0 

and 


0 . Then from 



= 0 . 


dUn , 

Sharp leading edge with — l jt 0.- At a sharp leading edge A = 0 
dc 

which reqiilres that (A)x=q “ ® from equation (15)* Hence (Z)^^.^ = 0 
and (0)x=o = ^ from equations (A3) and (Ak) . 

Substitution of the values of A and dA/dx from equation (I 5 ) 
into equation ( 29 ) and setting A = 0 result in 
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^ 1260 fir 2^ 

Wx=0 57 \Ti <3x 


(AS) 


The required expression for is then obtained from equation (A 3 ) 

V^/x=0 

after using equation (15), setting A = 0, and inserting. equation (a 8) 
which gives 




= iM (^^^)x=o 

^) y :=0 315 


(A 9 ) 


where for a .sbaip leading edge the quantity Rg^ would be designated as 

- 

ii.* V /x=0 

The initial value of B3_ is obtained as 




6300 

37 - 63O0aQ^ 


Oq + ^02(1 ~ 0) 


Ti* - Ti 

Tl* - \ 


(AlO) 


x=0 


by substituting the values of A and < 3 A/dx from equation (15) and cp/A 
from equation (27) into equation (30) and setting A = 0 . The required 
values of B2 and are evaluated from equations (A 5 ) and (a 6) with 

A = 0 . It is of interest to note that^ for q^ = 0 ^ B]^ = 0 from equa- 

tion (35) ^ and then equation (AlO) reduces to equation ( 4 o) with A = 0 . 


The initial valiie of ^ is obtained by rising the derivative of 

dX 

equation (15) in the derivative of equation (a 4 ) and setting A = 0 . The 
result is 
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after using equaticaas (a 8) and (AlO). The expression for is 

\<3x /x=o 

too lengthy to include here, particularly since this expression, although 
convenient, is not essential to the integration procedure. 


Sample Problem 

The present method is applied to the calculation of local heat- 

transfer rates on a circular cylinder at a stream Mach number of 6.9^ 

dT T„ 

= 0, p = 0° and 60°, and = 0.2 and 0.5. Typical results of 

this calculation are shown in figure 9 where Z, "X, 0, and B]_ are 
plotted against x. Note that in all cases A goes to 12 at some value of 
jE< - fx i= - corresponds to 90° back on the cylinder) as shown in fig- 

ure 9(a) . The calculation breaks down for X > 12 because of the incom- 
patible values of Z and — required to satisfy eq-uations (A3) and 

dx 

(Al) for X > 12. The physical reason for this is that the local chord- 
wise velocity within the boundary layer tends to exceed the external 

chordwlse velocity due to the large values of — ^ and — . A similar 

dx Tp 

phenomena has been observed in certain exact solutions of wedge-type 
flows (ref. 17)* 
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Figure 2.- Con^arlson of the calculated and approximate variation of F 

vlth 2V/T1 o = 0.70). 




MCA. IM L55F09 




O i 2 J 4 s 6 7 6 


Cross-How Mach /?umber , M^j = cos/} 

* ■ t^m ^ / 

Figure 3«- Con^arlson of experimental and theoretical values of velocity- 
gradient parameter at the stagnation line of a circular cylinder. 
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Figure 4.- Chordwise velocity profiles at tlie stagnation line as obtained 
from the present method and as given hy Crabtree (ref. 8). Both 
curves are for the conditions of zero heat transfer, Prandtl number 
uni ty, and spanwlse Mach number uniiy or = 1.2, and 

P1q/Pi* = 0.654. 
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Figure 5 .- Theoretical and experimental variation of heat transfer and 
recovery temperature with yaw angle on a circular cylinder at 
= 6 . 9 ^ and Ti* = 1200° R. 
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Figure J.- !Hie ratio of the calculated local heat transfer to the value 
at X = 0 on a circular cylinder at yaw angles of 0 ° and 6o° for 
= 6.9 and T;l* = 1200 ° R. 





Yah' anp/e , 

Figure 0.- [Qie theoretical and experimental variation of the average heat- 
transfer parameter with yaw angle on the front heif of a circular 
cylinder at 6.9 aiui « 1200*^ R. 
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(c) Variaticjn of Z and p with x at a yaw angle of oO . 

Figure 9*- Concluded. 
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